function xp = guanacos3(t, x)
global boxmuller i;
r= boxmuller(i);
a=2*(10^(-6));
s=0.01;
b=1*(10^(-6));
xp(1) = r*x(1) - a*x(2)*x(1);
xp(2) = -s*x(2) + b * x(1) * x(2);

xp = xp';